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Abstract. We have measured differential cross sections for pp elastic scattering with internal fiber targets in the 
recirculating beam of the proton synchrotron COSY. Measurements were made continuously during acceleration for 
projectile kinetic energies between 0.23 and 2.59 GeV in the angular range 30° < o .m. < 90° . Details of the 
apparatus and the data analysis are given and the resulting excitation functions and angular distributions presented. The 
precision of each data point is typically better than 4%, and a relative normalization uncertainty of only 2.5% within an 
excitation function has been reached. The impact on phase shift analysis as well as upper bounds on possible resonant 
contributions in lower partial waves are discussed. 



PACS. 25.40. Cm Elastic proton scattering - 
21.30.-x Nuclear forces 



13.75. Cs Nucleon-nucleon interactions - 13.85. Dz Elastic scattering 



1 Introduction 

Elastic nucleon-nucleon scattering is a process fundamental to 
understanding nuclear forces. Its precise experimental knowl- 
edge bears on such basic questions as the confinement of quarks, 
the limits to the validity of meson exchange models in the nu- 
clear regime and a basis upon which to test QCD inspired mod- 
els at intermediate energies. Its knowledge also forms the ba- 
sis of a broad range of applications in nuclear and heavy ion 
physics, e.g. as ingredients to models of reaction dynamics, ex- 
cited nuclear matter and transport phenomena. Consequently, 
many experimental and theoretical studies (see |1] and refer- 
ences therein) have been devoted to this subject. The data base 
has about doubled over the past decade Ol^QLu], an< ^ global 
phase shifts - a convenient tool of parameterizing existing ex- 
perimental knowledge - are on a sound basis up to about 1 GeV 
in kinetic beam energy 

Closer inspection shows, however, that the vast majority 
of the data is below 0.8 GeV of kinetic beam energy. Above 
that beam energy, the data base used to be increasingly sparse 
and somewhat uncertain in normalization, hampering both re- 
liable phase shift analysis and conclusions towards the physics 
of strong interaction. 
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The EDDA experiment was designed to provide pp elas- 
tic scattering data up to 2.5 GeV - precise and consistent in 
normalization - for all of the above purposes. In particular, it 
was meant to supply data for phase shift analysis up to that 
proj ectile e nerg y, and to test claims 0] and predictions, e.g. 
Ill (111 llll 2Ul 311. of dibaryonic resonances. To this end EDDA 
measured spin-averaged differential cross sections 1 14] as well 
as transverse analyzing powers lH5ll and spin correlation coef- 
ficients iTTrjll . utilizing the proton beams available at the cooler 
synchrotron COSY 1 17] at FZ Jiilich. 

An account of the results has already been given in fl4ll . 
based on roughly 40% of the data available now. The purpose 
of this paper is to describe experimental details and to present 
a more elaborate analysis. In particular, corrections for radia- 
tion damage of the target was not included in the analysis of 
Ref. [14] which lead to cross sections systematically smaller 
by « 3% for momenta above about 2 GeV/c. In addition, a re- 
fined separation of proton-carbon scattering events was accom- 
plished by pattern recognition techniques, resulting in reduced 
statistical errors. In total the statistical precision was improved 
by up to a factor of 3. 

The experimental setup is described in Sect.|2] while Sect. [3] 
focuses on a detailed description of the data analysis including 
event selection, vertex and angle reconstruction, background 
subtraction and normalization. Sect. |4] presents angular distri- 
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butions and excitation functions, Sect. [3]discusses their impact 
on phase shift analysis and hypothetical dibaryonic resonances. 



2 The EDDA experiment 

EDDA was conceived as an experiment using the internal beam 
of the cooler synchrotron COS Y IU7I1 . so as to measure contin- 
uous excitation functions during beam acceleration. It is placed 
at a beam waist in one of the straight sections of COS Y's race- 
track shaped lattice. With the protons recirculating with about 
1 .0-1 .5 MHz, inherently thin, polarized atomic beam targets for 
the measurement of spin observables can be utilized with suf- 
ficient luminosity. However, for cross section measurements as 
reported here, thin CH2-fiber targets have been used, as they 
allow for a continuous monitoring of the instantaneous lumi- 
nosity. 



2.1 Targets and Detector 

CH2-fibers were strung horizontally between the prongs of a 
fork, which could be moved vertically by a magnet-driven lin- 
ear actuator to put the fiber into (and out of) the COSY beam. 
The thickness of the fibers was chosen as a compromise be- 
tween long beam-lifetimes and sufficient sturdiness when ex- 
posed to the beam. A cross-section of 4 x 5 /im 2 CH2 proved 
to be a good choice, allowing a few 10 s protons recirculating 
in the COSY ring at about 1.5 MHz without target-failure. Tar- 
get fibers were prepared using a microtome to cut 4 /im fibers 
from a commercially available polypropylene (CH2) foil of 
5 /im thickness. To prevent build-up of charges due to ioniza- 
tion by the COSY-beam, the surface was made conducting by 
a thin (20 /ig/cm 2 ) aluminum coating, such that the target was 
grounded through the supporting fork made from aluminum. 

The background resulting from the carbon content (and the 
Al coating) needs to be subtracted offline; It was measured us- 
ing 5 /im thick carbon fibers suitably coated with Al. CH2- 
targets lose slowly hydrogen content upon beam exposure. There- 
fore, a reservoir of targets was kept inside the COSY vacuum 
to ensure essentially uninterrupted operation of EDDA during 
data taking. In first order, the relative normalization of the exci- 
tation functions is not affected by the slow hydrogen loss since 
the time period of about 2 s for measuring a complete excita- 
tion function during beam acceleration is small with respect to 
the hydrogen loss rate. 

The detector as used for cross section measurements is shown 
schematically in Fig.f^ A scintillator hodoscope surrounds the 
beam pipe downstream the target covering the angular range 
10° < 6 uh < 72° and subtending about 85 % of 4tt in the 
cm. for proton-proton elastic scattering. It consists of scintilla- 
tor bars and rings designed to measure the points, at which the 
scattered proton and its recoil partner traverse the hodoscope in 
terms of azimuthal and polar angle (ipi.2 and #i a b.i,2, respec- 
tively). The kinematic relationships between these angles (see 
Sect. 13. 6> are used to identify elastic scattering events. 

The inner hodoscope layer consists of 32 scintillator bars 
made of BC408 and read out on both ends via lucite light guides 
with Hamamatsu R1355 photomultipliers (PM) 1 18]. The scin- 
tillator bars have a triangular cross section (Fig. 13(a)). Their 
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Fig. 1. The EDDA detector (not to scale): target: fiber (CH 2 or C); 
B: scintillator bars; R: scintillator semi-rings; F: semi-rings made of 
scintillating fibers. 
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Fig. 2. Cross sectional view of adjacent scintillator bars (a) and 
semi-rings (b), crossed by two proton trajectories which can be dis- 
tinguished by the fractional light output of elements A and B. The 
first 9 semi-rings are made up of scintillating fibers as shown in the 
figure. 



overlap ensures that each charged particle originating at the 
target deposits energy in two adjacent bars when traversing the 
layer. The fractional light output from adjacent scintillator bars 
is used in the offline analysis (see Sect. 13. 51 to give the az- 
imuthal angle to a precision about five times better than would 
be possible on the basis of granularity alone 1 19]. 

The outer layer is composed of scintillator rings for angles 
$iab < 52°, split into left and right semi-rings to allow radial 
readout of the scintillation light in direction ±y. After total re- 
flection by 90° to the top or bottom of the detector, the light is 
collected via light-guides by Hamamatsu R1450 photomultipli- 
ers. The width of the rings varies along the beam (cf. Fig. 13(b)) 
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Fig. 3. Schematic setup of luminosity monitors. 

such that each proton trajectory crosses 2, in some cases 3 
rings. Each ring covers an interval A9 c . m , s» 5° in the center- 
of-mass frame for elastic pp scattering. Again the ratio of light 
outputs from adjacent rings is used to improve the resolution 
in (polar) angle over granularity. The method becomes less and 
less effective with increasing polar angle. Therefore, the scin- 
tillator rings were replaced by a double layer of (2 x 2 mm 2 
quadratic) scintillating fibers for angles 6\ a b > 52° as shown 
in the inset of Fig. 0(b) [20]. The latter cover about 10 cm of 
the inner layer close to the target, and they give - by their gran- 
ularity - about the same angular resolution as is obtained from 
the rings using fractional light output analysis. 



2.2 Luminosity Monitors 

Two independent, concurrent methods were applied to monitor 
luminosity as sketched in Fig. [5] Both monitors detect elec- 
trons emerging from the fiber target due to the beam-target 
interaction. The secondary electron monitor (SEM) consists 
of a fast amperemeter [21] measuring continuously the cur- 
rent of electrons from ground to target, replacing low-energy 
secondary electrons emanating from the target surface. High- 
energy (5-electrons from elastic proton-electron scattering were 
detected in two PIN diodes located at 6*i a b ~ 40° behind thin 
(250 /ig/cm 2 ) aluminum windows in small pockets of the beam 
pipe. Both rates scale differently with projectile momentum 
and must be corrected for this dependence when monitoring 
the luminosity consistently over the whole momentum range, 
see Sect.|^9] 



2.3 Cyclic operation 

Data are collected during synchrotron acceleration such that 
a complete excitation function is measured in each accelera- 
tion cycle. A typical cycle is shown in Fig. |4] After injection 
of some 10 7 protons with p = 275 MeV/c, acceleration starts 
with a momentum ramp of 1.15 GeV/c per second. At p w 
0.7 GeV/c, the target is moved vertically into the beam. The 
target remains in the beam until the maximum beam momen- 
tum of p ss 3.3-3.4 GeV/c is reached (flattop). Since the hor- 
izontal beam position is not constant during acceleration the 
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Fig. 4. Timing during an experimental cycle. Shown is the projectile 
momentum p (top), signal of the SEM luminosity monitor (middle) 
and trigger rate (bottom) as a function of cycle time. 



beam was steered across the fiber in the flattop (cf. Fig. [6), 
such that any effect relating to the position along the fiber tar- 
get could be studied offline. During deceleration the beam is 
lost and the synchrotron prepared for the next cycle. The SEM 
signal indicates that the luminosity roughly remains constant 
during beam acceleration since beam losses are compensated 
by the increased beam-current when the revolution-frequency 
raises from 0.5 to 1.5 MHz . It should be noted that the beam 
emittance growth due to small angle scattering in the fiber tar- 
get was nearly co mpen sated by the adiabatic damping due to 
beam acceleration [22]. During the flattop the luminosity de- 
creases with a lifetime in the order of 5-10 s. 

After retraction of the target, a laser monitor system (cf. 
next section) is activated for 2 s in the flatbottom part of the 
cycle. For calibration and beam diagnostics, alternative cycle 
modes are applied, e.g. measurements at fixed momenta with 
longer flattop for detector calibration, cf. Sect. 13.51 



2.4 Laser Monitor System 

The EDDA-detector was monitored continuously by a laser 
monitor system. Light pulses are generated at a rate of 20 Hz 
by a nitrogen laser which drives a dye laser. This produces light 
pulses with mean wavelength (A = 425 nm) and pulse shape 
characteristics similar to those generated in the scintillators by 
real particles. They are fed through individual fibers into the 
light guides of all scintillator elements to illuminate all photo- 
tubes simultaneously for an online detector control. Pulse-to- 
pulse variations in laser intensity are accounted for by normal- 
ization to photodiodes receiving a constant fraction of the dye 
laser light. In offline analysis the laser events are used to de- 
tect long term gain drifts or jumps between detector calibration 
runs and to deduce necessary corrections. 
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For very low energy ejectiles the produced light may ex- 
ceed the linear range of the PM-response. To this end, separate 
laser runs are performed with the pulse intensity being stepwise 
decremented by a set of calibrated optical filters. By compar- 
ison to the photodiode reference the functional dependence of 
the digitized collected charge from the PM on the light inten- 
sity is mapped out. Offline these correction functions can be 
used to linearize the response of all photo-tubes individually 
(EH and Sect.|33}. 

2.5 Trigger and Data Acquisition (DAQ) 

The structure and granularity of the scintillator hodoscope re- 
flect the signature of pp elastic scattering events, namely (i) 
coplanarity 

1^2-^1 = 180° (1) 
and (ii) kinematic correlation of scattering and recoil angle, viz 

tan6> lab ,i • tari6»i abi2 = -j— , (2) 

Tc.m. 

where 7 c . m . = yl + T p /2m p c 2 denotes the Lorentz factor of 
the pp center-of-mass motion as a function of beam kinetic en- 
ergy T p . Both conditions are used to define a fast online trigger. 
The coplanarity, Eq. Q, is verified by requiring a coincidence 
of two scintillator bars just opposite to each other with respect 
to the beam. Light attenuation along the bars is taken care of 
by adding the readouts from either end. The different flight and 
light transport times as well as electronic constraints require a 
sufficiently wide (70 ns) coincidence gate. Monte Carlo sim- 
ulations confirm that the width of the bars guarantees a trig- 
ger even for deviations of the vertex position by as much as 
±1.5 cm from the nominal COSY beam axis. 

The trigger condition for elasticity, Eq. correlates each 
semi-ring on one (e.g. left) side of the beam with one or more 
on the opposite (e.g. right) side. These correlations, however, 
depend on projectile energy T p as well as the transverse pro- 
file of the beam-target overlap. If the range of rings allowed in 
coincidence is sufficiently wide to be valid for the whole mo- 
mentum range covered during beam acceleration only a poor 
suppression of quasi-elastic scattering in the carbon-nuclei is 
achieved. For this reason programmable logic modules have 
been used to implement the coincidence, such that they can be 
reprogrammed within 5 ms while the beam is accelerated. A 
gain in trigger efficiency by a factor of two (cf. Fig. [5jl was 
achieved by reprogramming three times during acceleration, 
sacrificing 1 % of data acquisition time. The times when the re- 
programming occurred was changed frequently to avoid empty 
spots in the final excitation functions. 

Data are acquired whenever the target is in the beam or the 
laser system is triggered. Data are processed with conventional 
CAMAC modules for timing, pulse-height and logic signals in 
conjunction with a VME based event-builder. An event is writ- 
ten on tape when Eqs. Q and (|2ji are fulfilled by at least one 
pair of scintillator bars and semi-rings, respectively. In addi- 
tion, each event contains the cycle time from a COSY clock 
controlling the accelerator operation, the collected charge from 
all scintillator elements and the timing signals from both ends 
of the scintillator bars. 
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Fig. 5. Scatterplot of events taken with a carbon target during accelera- 
tion. By threefold reprogramming (dotted horizontal lines) of the kine- 
matic trigger, the reduced acceptance (solid lines) suppresses uncorre- 
cted background that cannot be avoided with a fixed setting (dashed 
lines). 

There are three additional trigger sources, namely for PIN 
and SEM luminosity monitors and for the laser system, which 
are merged with the main trigger. All four independent trigger 
sources are counted, together with a 10 MHz clock and pulses 
counting the instantaneous frequency ^rf of the COSY cavi- 
ties, in scalers which are read out in time increments of 2.5 ms. 
Comparing the number of stored events with the scaler of the 
corresponding trigger source yields the dead time fraction r(p) 
of the data acquisition system as a function of beam momentum 
(i.e. cycle time) that can be corrected for in beam momentum 
intervals as narrow as 3 MeV/c (cf. Sect. 13.1 IV Our measure- 
ments were performed with time-averaged luminosities of typ- 
ically 5 • 10 29 cm~ 2 s _1 and singles rates of up to 40 kHz. 



3 Offline analysis 
3.1 Data Samples 

The data were collected in two production runs of two and 
three weeks, separated by half a year, with a yield of about 
15 TO 6 and 22 10 6 elastic events, respectively. After publica- 
tion of the results based on the first production run 0, small 
modifications concerning e.g. the treatment of artefacts due to 
the detector granularity have been implemented in the analy- 
sis, advanced methods to reduce non-elastic background were 
developed, and a correction for radiation damage of the CH2- 
target introduced. Furthermore, the multitude of different run- 
ning conditions, have allowed additional consistency checks 
and both data sets have been (re-)analyzed for the present work. 

First, the data have been divided in 17 so-called data sam- 
ples, which have been analyzed separately. Each sample com- 
prises typically 1.5 . . . 4 • 10 6 elastic events and consists of two 
data sets, one taken with a CH2- and a C-target each. The rela- 
tive luminosity of the CH2- and C-target data, measured inter- 
leaved by changing targets every few hours, was chosen online 
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to yield about the same statistics for proton-carbon scattering 
events. Each data sample contains data 

1. from only one CH2- and one C-target, 

2. with exactly the same trigger setting (e.g. reprogramming 
times) 

3. and the same detector and COSY-beam setup. 

By treating the individual results for these samples as separate 
measurements many consistency checks could be made prior to 
combining them to the final results (see Sect. 13.121 . 



3.2 Overview 

For N pp detected elastic scattering events within a solid angle 
bin A9 c . m . and a momentum bin Ap, centered at (p, c . m .), the 
differential cross section is given by 



da 



(p,i 



N PP (p 



pp\Fi u cm. 



AS2.n{p,o c . m )-L B {p).{\-T{p)y 



(3) 



Here, Lh denotes the absolute luminosity and r the DAQ dead 
time fraction for the respective momentum bin. The detection 
efficiency 77 accounts for the elastic scattering events discarded 
by the trigger or offline cuts, mainly due to secondary reactions 
of scattered protons in the experimental setup (Sect. l3~71 . 

For elastic proton-proton scattering we detect two protons 
and a rough selection of events is based on detector multipli- 
cities. To this end, we subdivide the data in event classes based 
on the number of hits in the scintillator elements of the ho- 
doscope and keep only those of low multiplicity. The scatter- 
ing angle 6 C . m . is reconstructed from the vertex position and 
the interception points of the two trajectories, determined from 
the energies deposited in the scintillator elements. The vertex 
position moves along the target fiber during beam acceleration. 
Applying a kinematic fit, it can eventually be determined event- 
wise. The corresponding concepts are described in Sect. 13.51 
and were applied to measurements with both types of target 
fibers. 

For background reduction, a combination of an event-wise 
selection based on Eqs. Q and (0 and a statistical subtraction 
of the carbon fiber data, normalized to the same luminosity, is 
accomplished in a way that - in first order - inelastic pp inter- 
actions are corrected for as well. Details are given in Sect. 13.61 
For all excitation functions, the relative luminosity is provided 
by the two luminosity monitors. The determination of the ab- 
solute luminosity L(p) in this experiment is not as accurate as 
in some high-precision external experiments, so that we nor- 
malize all excitation functions at one single momentum p to 
published data (cf. Sect.l3~9l. 



1 . Calculating the efficiency r\ for accepting elastic pp-scatter- 
ing events, this includes effects of the trigger as well of all 
software cuts used in the analysis. 

2. For modeling the contribution of inelastic pp-reactions in 
order to estimate the background contribution to accepted 
pp-elastic events. 

3. To calculate the detection efficiency of the PIN-diodes for 
(5-electrons used for measuring the relative luminosity. 



For electromagnetic interactions we use EGS4 |36] and for 
hadronic reactions MICRES |25] or phase-space distributions 
as event generators. In MICRES a number of event generators 
are used: To model elastic scattering cross sections distribu- 
tions according to recent phase-shift analysis |7] are generated. 
Inelastic scattering of hadrons on nucleons and nuclei for ener- 
gies up to 5 GeV are described by tracking the impinging nu- 
cleon through the nucleus. Secondary particles including pions 
and hadrons are traced and their interactions with nucleons are 
described with known, interpolated or estimated cross sections. 
Momentum, energy and charge are fully conserved. In this in- 
tranuclear cascade the nucleon-nucleon interactions are either 
quasi-elastic or inelastic, which is modeled by exciting reso- 
nances, which may subsequently decay and eventually produce 
mesons. All nuclear resonances contributing more than 2 % to 
the total inelastic cross section are taken into account. Details 
of this event generator MICRES are given elsewhere [25]. 



3.4 Multiplicity Cuts 

For a first classification of an event, a group of n hits in adjacent 
scintillator bars or semi-rings is called a cluster of size n. Due 
to the geometrical overlaps (cf. Fig.|2) cluster sizes n = 2 and 3 
dominate. An event can then be characterized by its cluster pat- 
tern (Nb, Nl, Nr), where is the multiplicity (i.e. the num- 
ber) of clusters in the layer of scintillator bars, and (N/{) 
denotes the respective multiplicity for the left (right) group of 
semi-rings. An ideal elastic pp event reveals the pattern (211). 
The majority of elastic scattering events is of this type. There 
are, however, several patterns of higher cluster multiplicity that 
originate from elastic scattering due to e.g. secondary reactions 
in the scintillator hodoscope, accompanying ^-electrons or ran- 
dom coincidences. By comparing the CH^and C-data sets the 
patterns containing non-negligible numbers of elastic events 
could be identified and we kept eight patterns, namely (211), 
(311), (411), (221), (212), (321), (322) and (222). For a projec- 
tile energy T p = 1.5 GeV, these patterns exhaust 98.4% of all 
triggered events, with the ideal (211) pattern already contribut- 
ing 89.2 %. All other patterns are discarded. 



3.3 Monte Carlo simulations 

Experiment and analysis are accompanied by a comprehensive 
Monte-Carlo simulation including the details of detector ge- 
ometry and materials. Particles are tracked through the detec- 
tor, taking into account energy and angular straggling as well as 
hadronic and electromagnetic secondary interactions. Results 
of Monte-Carlo studies enter mainly threefold into the analy- 
sis: 



3.5 Trajectory and vertex reconstruction 

EDDA primarily detects the position of the protons on cylin- 
ders surrounding the beam pipe. The coordinate along the cir- 
cumference R ■ 4> is determined from the scintillator bars (with 
R = 16.4 cm) and the z position along the beam by the semi- 
rings. The resolution can be enhanced considerably beyond that 
given by the detector granularity by evaluating the charge-inte- 
grated photomultiplier-signal S of all scintillator elements, usu- 
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ally two, hit by the particle. This can be related to the path- 
length within the scintillator element Ax and to the average 
specific energy loss along this pathlength (dE/dx) by 



' dE\ 

S(6 lah ,<t>)=G- [ —J -Ax 



(4) 



For most projectile energies and scattering angles 0\ a b covered 
in the experiment, the energy loss can be considered constant 
along the whole path through adjacent, overlapping scintilla- 
tor elements. In this simple case the gain factor G describes 
the conversion of deposited energy to the electronic signal and 
encompasses (i) light transportation and attenuation in the scin- 
tillator and light guide, and (ii) photomultiplier (PM) response. 
The former effect introduces a dependence on the particle's 
point of incidence G(#i a b, <j>) which is mapped out experimen- 
tally by pp elastic scattering data taken at fixed momentum 
p = 2.7 GeV/c, where the detector acceptance is at its max- 
imum. 

For the latter, a laser-light (cf. Sect. 12. 4> fed into the light 
guides was used to deduce corrections for long-term gain drifts 
and nonlinearities close to the upper limit of the PM's dynamic 
range. 

The spatial resolution in <fi (9\ a h) can then be enhanced be- 
yond the granularity given by the number of scintillator bars 
(semi-rings), by interpolating the point of incidence using the 
pulse heights S%, S2 of neighboring scintillators traversed by 
the charged particle. The ratio 



Qi 



S2/G2 — S1/G1 
S2 1 G2 + S\jG\ 



(5) 



varies according to Eq. monotonically with Ax2 — Axi 
between +1 and —1 along the interval of overlap, see Fig. [2] 
Details of the algorithm are given in 1 19]. The resolution ob- 
tained for an ejectile is improved by a factor ps 5 to 6(f) « 1.9° 
(FWHM) and 58 c . m . » 1.0° (FWHM). 

For ejectile energies below 400 MeV two corrections are 
important: First, the specific energy loss dE/dx can no longer 
be considered constant along the path of a particle. Assuming 
that the particle is an elastically scattered proton, so that the 
kinetic energy is known and taking into account the stopping 
power of all material in front of the detector element under con- 
sideration a correction can be deduced. Secondly, the amount 
of produced light is affected by nonlinearities in PM response. 
Therefore, it is linearized using information obtained with the 
laser monitor system (cf. Sect. 12. 4> . With these corrections a 
good angular resolution is maintained except for the lowest en- 
ergies, where systematic deviations of the reconstructed posi- 
tion arises which will be treated as a systematic error. 

In order to obtain the scattering angle the reaction vertex 
must be known. The vertical and longitudinal position are fixed 
by the fiber position and are thus the same (within a few ^m) 
for all events. The horizontal vertex position is smeared out 
by the width of the COSY beam and should be determined for 
each event separately. As a first step the first moments of the 
vertex distribution, e.g. the mean values (x v ,y v ,z v ) and the 
horizontal width 5x v are deduced from the data, by exploiting 
the coplanarity of elastically scattered protons with the beam 
or the kinematic correlation. 
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Fig. 6. Variation of horizontal vertex parameters and vertical beam 
position y v as a function of cycle time. The vertical line shows the 
begin of the flattop where the beam is swept on purpose horizontally 
over the fiber target. x v is the beam position, Sx v the beam width 
folded with the detector resolution, x' b the horizontal beam angle and 
y v the vertical target position. The first two points in y v show the 
settling of the target after it has been moved in. 



The vertex position x v , y v in the plane perpendicular to the 
symmetry axis of the detector can be determined from the pairs 
of points where the prongs intercept the scintillator bars. Due 
to the coplanarity of two body interactions the projection of 
the line connecting the two points into this plane must cross 
the vertex. The elastic scattering events are evenly distributed 
in cf>; the coordinates x v , y v are obtained as parameters from 
a fit to all events for a given projectile momentum interval 
Ap (details are given in flfx ). The method can be extended 
to projectile beams tilted against the detector axis. By compar- 
ing events with equal (6> c . m . ps 90°) or different (0 c .m. < 90°) 
lab scattering angles yields the directions x' b and y' b as well (cf. 
Fig-EJl. The validity of this method was tested and verified with 
a beam of fixed momentum that was steered in a controlled way 
horizontally (x) along the target fiber (3.6 s < t < 4.8 s in 
Fig- EJl. The horizontal vertex position and the tilting angles of 
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the beam turned out to vary during acceleration to a non negli- 
gible extent. These parameters, however, showed a remarkable 
long term stability over many days. Finally, the vertex position 
z v in beam direction is calculated from the interception coor- 
dinates z\, Z2 in the ring-layer by making use of the kinematic 
relation Eq. (0 with 7 c . m . being known from the instantaneous 
projectile momentum and tan#i a b,i = R/zi, where R is the 
mean radius of the semi-ring layer. This way, the vertex coor- 
dinates x v , y v , z v are determined with an accuracy ± 0.5 mm; 
the width Sx v « 3 mm (FWHM) reflects the horizontal beam 
width folded with the experimental resolution. 

With the mean vertex position and the position of the hit 
in the detector known, the polar and azimuthal scattering an- 
gles are calculated, taking into account the small correction 
due to the tilted beam. With the reconstruction algorithm de- 
scribed above, however, the horizontal vertex position cannot 
be determined eventwise. To improve the resolution of the final 
cm. polar scattering angle 6* c m , it is taken from a kinematic 
fit, constraint by elastic scattering kinematics (i.e Eqs. and 
(j2ji). In the fit x v is treated as a free parameter; its distribution 
is compatible with the previous result for 5x v . 

For events with higher multiplicities, i.e. pattern other than 
211, the ambiguities when correlating hits from the bar and the 
ring layer can be resolved almost entirely by the time-difference 
of the signal read out at the down- and upstream sides of the 
scintillator bars, which yields a position information along the 
bar of about 5 cm (FWHM) and must match the position of 
the struck rings. If ambiguities remain, the combination of hits 
which best matches elastic scattering kinematics, as outlined in 
the next section, is selected. 



3.6 Event Selection and Background Subtraction 

For each event we have at least two (i = 1,2) reconstructed 
tracks with angles (6\ a b,i,4>i) mat ^ now subject to additional 
cuts to separate elastic pp scattering events from inelastic reac- 
tions or background originating from proton-carbon reactions. 
For this purpose one first calculates the cm. scattering angle 
using pp elastic scattering kinematics: 



1 - cos( a) 
0.01 0.05 0.1 



kin. deficit a (deg) 
5 10 



2 arctan (7 c . m . tan#i ab ,i) 



(6) 



Note that <fi — <p\ a \> — 4> c . m .. For elastic events both particles 
should be emitted back-to-back in the center-of-mass system. 
The angular deviation a from this perfect 180° correlation, 
subsequently called the kinematic deficit, is shown in Fig.0for 
CH2- and C-target data at two different beam momenta. The 
kinematic deficit is used twofold: First, for the few events with 
more than two prongs the combinatorial ambiguity is removed 
by selecting the combination with minimum a. The resulting 
probability distributions N(j>, a) show that elastic pp events 
clearly stand out on a smooth, monotonic background. A com- 
parison to the distribution obtained with a pure carbon fiber 
target under identical experimental conditions shows that this 
background can essentially be associated with quasi-free scat- 
tering on bound protons and inelastic p-C interactions. A dis- 
play of events more closely related to the conditions of copla- 
narity, Eq. Q, and elasticity, Eq. (|2}, is shown in Fig.[8]where 
the elastic peak extends about two orders of magnitude above 



10 rv 



p= 1.45 GeV/c 



p = 3.00GeV/c 



p= 1.45GeV/c 



p = 3.00 GeV/c 



10 20 30 0.005 0.01 0.015 0.02 



kin. deficit a (deg) 



1 - cos( a) 



Fig. 7. Distribution of the kinematic deficit a for two beam momenta. 
Shown are events from CH2(closed) and C-targets (open symbols). 
The elastic pp events are indicated by the shaded area. The upper limit 
a cut (cut) for accepted elastic scattering events and the range (norm) 
used for normalizing the C-target data to the CH2-data are shown. The 
solid angle associated with a rises with cos a, so that the number of 
events at very small a drops. The figures on the right show that indeed 
the elastic peak is at cos a = 1 as expected. Note the logarithmic 
scale. 



the background. However, the combination of those two re- 
quirements in one variable, a, is more effective for quantitative 
background treatment. 

In a first step, a projectile momentum dependent cut a < 
CKcut(p) is applied in a way that - based on Monte-Carlo simu- 
lations and comparisons to the C-target data - almost no elas- 
tic events are discarded. The remaining events are sorted into 
two-dimensional arrays iVcH 3 (P>0c.m.) an d Nc(p,Q c .m.) with 
the option to adjust the bin- widths Ap and A6 c . m . to statisti- 
cal requirements. 



3.6.1 Statistical Background Subtraction 

Background subtraction is performed statistically: A factor C(p) 
is deduced from the normalization in the interval = a cu t + 
4° < a < a cu t + 11° = a max of the a-distribution 
NcH 2 (p,ct) of a CH2-target data set to that Nc(p,ct) of the as- 
sociated C-target data set, i.e. 



£(p) 



a m i n <a<a n 
^0(P)la tata <a<a n . 



(7) 



where events have been integrated over all scattering angles. C 
can be viewed as the luminosity with respect to proton-carbon 
scattering of the CH2-target relative to the C-target data set. Its 
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Fig. 8. Distribution of angle integrated events obtained with a CH2 
(top) and carbon (bottom) fiber target at p = 2.25 GeV/c. Note the 
logarithmic scale. 



momentum dependence is caused by the different COSY beam 
lifetimes for the two targets due to their different thickness. C 
was determined for all contributing multiplicity patterns sepa- 
rately, to account for slightly modified contributions from ac- 
cidental coincidences caused by the luminosity difference. The 
background can now be subtracted statistically viz 



N w (p, 9 c . m .) = N CK2 (p, 9 c . m .) - £(p) • N c (p, 6 c . m .) (8) 




p = 3300 MeV/c 



a pp elastic 
□ pp inelastic 
■ pC 



5 10 15 20 25 30 

kinematic deficit a (deg) 

Fig. 9. Simulated distributions of the kinematic deficit a for p-C re- 
actions (solid dots), elastic (triangles) and inelastic (squares) pp scat- 
tering. The vertical line indicates the position of a cut . 



3.6.2 Contributions of Inelastic pp-Reactions 

Elastic scattering exhausts at least half the pp total cross sec- 
tion in the COSY energy range. Except for pp-^d7r + all inelas- 
tic reactions have at least 3 particles in the final state and there- 
fore the probability to look like an elastic event is expected 
to be small. Since these reactions cannot be measured sepa- 
rately by our experiment we rely on Monte-Carlo (MC) sim- 
ulations, using the available information on the respective to- 
tal cross sections. Unfortunately, phase-space distributions are 
practically unknown so that models have to be used in MC 
event generation. For this purpose distributions either accord- 
ing to phase-space or generated by MICRES, where inelastic 
reactions proceed through excitation of resonances, are used. 

Simulations show that the main inelastic contribution orig- 
inates from the channels pp^ pn7r + , — > pp7r + 7r~ and — >pp7r°. 
The cuts on the cluster pattern and the correlation mismatch 
a reduce this background to less than 1.5% of the (angle inte- 
grated) elastic yield, whereas the corresponding number for the 
p-C background is up to one order of magnitude higher; both 
increase with beam momentum p. 

Their simulated a-distributions are shown in Fig.[9]for the 
worst case, i.e. at high momenta. Both background distribu- 
tions are similar in shape. Therefore, the statistical subtraction 
with its normalization performed in the angular range a > a cu t 
takes care of the main part of pp inelastic events with a < a cu t, 
too. A closer look, however, reveals differences in shape that 
will give rise to a systematic background error which has been 
studied in more detail. 

For this purpose, Ncn 2 (P> #c.m.) has been numerically com- 
posed from three contributions, namely 



There are two concerns when applying this method: first, the 
CH2-targets are aluminum-coated and it must be verified that 
p-C and p-Al reactions in this region of phase-space are indis- 
tinguishable. This was proven by taking data sets with plain and 
aluminum coated C-targets. The carbon-subtracted pp elastic 
scattering yields is the same within statistical errors. Applying 
the statistical subtraction scheme to C+Al- and C-target data 
sets (instead of CH 2 and C) produced results for "N" com- 
patible with zero. A second concern is that inelastic proton- 
proton scattering populates the a-distribution in the region where 
C is determined and will therefore systematically distort the re- 
sults. To shed light on this issue Monte-Carlo simulation were 
performed to mimic the effect of inelastic pp reactions. 



^ch 2 - ATpp + N£ + < a) (9) 

where the elastic contribution N° l p is based on our experimen- 
tal results 1 14]. The second term is Monte Carlo generated and 
the third is an experimental data sample Nq obtained with a 
carbon target, properly normalized to the observed ratio of pp 
to p-C events. This set is then corrected for background with the 
statistical method, Eq. (|8j, by using a second, statistically in- 
dependent carbon measurement Nq . From the result the term 
A^pp used before in Eq. (|9} is subtracted. The remaining dif- 
ference D(p : 9 c . m ) reflects the incorrect treatment of the in- 
elastic pp component A/I"", and turns out to be systematically 
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40 60 



Fig. 10. (a) Ratio of events from carbon- and CBh-targets of one data 
sample after all cuts have been applied and the carbon data have been 
properly normalized. The solid line shows the acceptance cut on the 
cm. scattering angle. In (b) the same data are shown when in addition 
event selection using fuzzy-logic methods is utilized. 



positive, indicating that the inelastic pp contribution was in- 
sufficiently subtracted as is expected from Fig. [9] The fraction 
|D/iVpp| is negligible (i.e. less than 1%) except for the large- 
angle (9cm. > 70°) and high-momentum (p > 2.5 GeV/c) 
region where the contribution may reach about 5.5%. Because 
only the total cross sections but not the angular distributions 
of most inelastic pp reactions are sufficiently well known, no 
attempts have been made to correct this deviation and the max- 
imum effect obtained with different models is taken as a sys- 
tematic uncertainty. 



3.6.3 Hydrogen Absorbed in Carbon-Fibers 

A close look at the distribution for events measured with carbon- 
fiber targets reveals a tiny peak at the locus of elastic pp scatter- 
ing, cf. Fig. [8] which could be due to traces of hydrogen or wa- 
ter absorbed in the carbon fiber. Indeed does this peak show an 
angular distribution in agreement with the expectation for elas- 
tic pp scattering. The background subtraction following Eq. (|8j 
will then reduce the number of elastic pp events of the P-CH2 
measurement systematically wrong. A careful study of this re- 
duction showed that it is for all targets smaller than 0.5% and 
does not vary with projectile momentum. Therefore, it is equiv- 
alent to a reduction of the effective CH2 target thickness and its 
effect cancels completely when finally all excitation functions 
are normalized at one momentum to some precise reference 
cross-section, cf. Sect. 13.9.41 



3.6.4 Statistical Errors 



The statistical error associated with the background subtrac 
tion is roughly given by SN pp = \ l N, 



CH 2 



C 2 N C where the 



error of C, which is a minor contribution, has been left out for 
the sake of simplicity. Thus any contribution from scattering 
off carbon that can be eliminated by cuts prior the subtraction 
procedure (Eq. (|8jl) will reduce SN pp . In Fig. [K)](a) the frac- 
tion of events in N cli2 attributed to the carbon-content of the 
target is shown. It rises with momentum and angle from 5% 



up to 25% and unfortunately is largest where scattering rates 
are the lowest. Therefore, it was investigated how the statistical 
method could be complemented by an eventwise classification 
scheme to reduce N c prior to being subtracted. To this end 
information on each event besides the measured angles, like 
energy losses, can be used. When testing methods commonly 
employed in pattern-recognition, like self-organizing maps or 
artificial neural networks, we found an approach based on so- 
called fuzzy-logic to be particularly helpful. It allowed a large 
reduction of the number of proton-carbon scattering events to 
be subtracted statistically as shown in Fig.^|(b), and will be 
discussed in detail in the next section. 



3.6.5 Background subtraction using Fuzzy-Logic 

So far only the information on the proton angles have been used 
to identify elastic scattering events. In addition, the EDDA de- 
tector provides energy-loss information in the bars and semi- 
rings and timing information from the bar layer, which can be 
put into use to test their compatibility with elastic pp-scattering 
events. However, quasi-elastic scattering from nucleons in the 
carbon-nucleus as well as inelastic pp reactions are not cleanly 
separated in these variables, so that simple cuts are insufficient. 
Here, methods developed within the framework of so-called 
"fuzzy-logic" [27, 28] have been shown to be very powerful 

& 

For each event, described by a vector q of n measured or 
reconstructed quantities qi (i = 1 . . . n), we must construct a 
decision-function d(q) to decide if it is accepted (d(q) = 1) or 
not (d(q) = 0) as a candidate for an elastic scattering event. 
When using simple cuts on one or two-dimensional projections 
of q applied in sequence most of the correlations between the 
observables qi are lost, since one variable outside one cut is 
sufficient to discard an event. In the fuzzy-logic approach we 
define for each observable qi a so-called member-function Pi : 
qi h—s- [0,1], which maps the parameter space of qi onto the 
interval between zero and one. It is designed such, that events 
with a high (low) probability of being elastic yield values pi 
close to 1 (0). For this purpose the member function must not 
be a probability in the strict sense and we use 



PMi) = 1 - 



(10) 



To obtain Pi with sufficient statistical precision, the data are 
binned in qi. N c and A r CHa are the number of events within 
a certain bin centered at qi for this quantity, whereas all other 
quantities qj with i ^ j may have arbitrary values. Note, that 
here and in the following it will be assumed that the CH2- and 
C-target data samples have been properly normalized to the 
same luminosity. 

These member-functions are then combined to a combined 
"probability" by the 7-operator 127T.I28T1 



p(s)( ( 



1-3 



(id 



mapping the parameters space in q onto the interval [0, 1]. It 
features a, yet to be fixed, parameter g £ [0, 1], to the effect 
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Fig. 11. Distribution of events in the (P^},Pj^)) plane for a CH2 (a) and C-target (b) normalized to the same luminosity. In (c) the difference 
attributed to pp-scattering is shown, and in (d) the "probability" P e \, as defined in Eq. i 1 31 . The dashed line shows a typical contour for accepted 
events based on the condition P c \, > P cut which corresponds to a loss of 3% of elastic events (e) while reducing the background by about 
50 %. The spectra shown correspond to events around p w 2.1 GeV/c and are integrated over all scattering angles. 



that the 7-operator specializes to the well-known logical AND 
(g = 0) or OR (g — 1) operators in the limiting cases. By this 
procedure we have mapped the information contained in q to a 
single variable P^> between and 1. 

In this experiment we used n=5 experimental quantities qc 

a: 
A<f>: 



A 



Tz 



The kinematic deficit as defined in Sect. 13.61 
The difference of the azimuthal angles of the two 
tracks, 

The correlation of the time-of-flight (TOF) T x - T 2 
and the z -position difference in the bar layer of the 
two charged tracks: 



An 



(Ti - T 2 ) - (3 c .m. (zi - z 2 ) , 



(12) 



where /3 c . m . — sJT v j (2m p + T p ). This quantity is 
close to zero for elastically scattered protons as re- 
quired by kinematics. 
AEi(2)'- The energy loss in the bar layer for the particle scat- 
tered to the left (right) of the beam. 

These were combined to two vectors 5( a )= ( a > A(f>, At z ) and 
q(b)= (AEi,AE 2 )- This allows us to map each event onto a 

unit square spanned by pijjv (<7( a )) and Pxj) (<?(;,)) as shown in 
Fig.[^for g = 0.5. The elastic events clearly stand out in the 



upper-right corner in the (P$ ,P$) plane, whereas reactions 
on carbon - considered as background - are spread out over a 
larger region. Increasing the value of g spreads out the back- 
ground over a larger area, whereas a decrease does the same to 
the elastic scattering events, so that g = 0.5 appeared to be a 
good compromise. 

In principle one could now define the decision function 

d(q( a ) 1 Q(b) ) by drawing a two-dimensional contour on the pP? 



vs. P$ distribution. However, the signal (pp elastic scatter- 
ing) and the background (p-C reactions) do overlap, so that any 
reduction of the latter will always be accompanied by some 
unavoidable loss of signal. However, in large regions of the 
(PW,P$ )-plane the losses will be very small. This can be 
quantified by defining P c \ ., again viz 



P eh (P^P^)^l- 



(9h 



(b) 



>(9)> 



(13) 



(b) 



If we sort the data in 20x20 bins in the (PW,P$) plane (cf. 
Fig. \^\ (d)) we can view P e i. as a probability that an event 
falling into a certain bin is from elastic scattering. Then, it is 
optimal to discard events with a low P e i., i.e. defining a deci- 
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Fig. 12. Spectra of the kinematic deficit a for two beam momenta 
ranges (±150 MeV/c) with (dashed line) and without (solid line) ap- 
plication of background rejection using fuzzy-logic. Shown are the 
data from CH2- (left) and C-targets (center, normalized to the same 
luminosity) as well as the difference (right) attributed to the hydrogen 
contents. Note that a cut at 7.5 ... 8° (top) or 4.5 ... 5° (bottom) on 
a has already been applied. 



sion function: 



d(g) 



1 for P e i. (q) > P cut 
for Pel. (9) < ^cut 



(14) 



where P C ut(p,S c . m .) is to be chosen momentum and angle de- 
pendent. Since the background originating from the carbon has 
been measured separately, we can now quantify both the back- 
ground reduction as well as the loss of elastic events as a func- 
tion of P cu t- In Fig.^2( e ) we show a typical example. How- 
ever, this procedure is done for every bin in (p, # c . m .) and ev- 
ery data sample individually. One obtains a sizeable reduction 
of background events (C) for minor losses of elastic scattering 
events (H). We decided to sacrifice 3% of the elastic scatter- 
ing events, allowing for a reduction of the background of about 
60%. A typical contour of accepted events is shown as a dashed 

line in Fig. [^(d). The events with P$ are populated 
mainly by events where the time information of the scintillator 
bars is ambiguous due to multiple hits. This occurs for some 
elastic events as well, so that this class had to be kept. 

In order to determine P C ut(p,#c.m.) with high statistical pre- 
cision we determined its value for 0.3 GcV^ wide momentum 
and 10° wide # c m -bins. Its value is typically around 0.5 and 
smoothly varying with the beam momentum and cm. scatter- 
ing angle. When applying the cut to the final data, its value has 
been interpolated for the correct momentum and angle. 

To summarize: by transferring the detector information - 
by methods adopted from fuzzy-logic - to two variables with 
values between and 1, any event can be mapped onto the 
unit square. Since the data are binned we group the data in 

400 classes according to their P^'and P$ values. Finally we 
select a certain subset of these classes by requiring that the re- 
spective P c i. is above a certain threshold P cu t, chosen to result 
in only a small, well-defined loss of elastic events. This reduces 



Fig. 13. Relative reduction of the statistical error from carbon sub- 
traction when applying the cut from the fuzzy-filter approach (dashed 
line) and the reduction on the total uncertainty (solid line) when all 
other sources of errors are included. 



the number of events, N CH (p, c . m .) and N c (p, 6* c m ), enter- 
ing Eq. to correct for the contribution arising from carbon- 
contents of the target (Sect. 13.6. Q . Since less events are sub- 
tracted (Figs. [|2l(b) and !12> the corresponding statistical un- 
certainty is considerably reduced (cf. Fig.1131. most notably at 
high momenta and large cm. scattering angles where elastic 
cross sections are smallest. 

Since the resulting loss of elastic events is independent of 
beam momentum and angle and its absolute value can be de- 
termined with 0.5% precision (e.g. 3% ± 0.5%) this is not off- 
set by a large increase of the systematic uncertainty. To this 
end a comparison of the final cross sections with and without 
pre-selection of events using fuzzy-logic (Fig. I14> shows no 
systematic differences. The precise value 0.5 of the parame- 
ter g introduced in Eq. (II Q is not important, however, values 
close to or 1 are to be avoided, since either the elastic events 
are concentrated in only a few classes or spread out over all 
classes. Then P e j, can not be determined with sufficient preci- 
sion for each class and the net effect on background reduction 
is diminished. 



3.7 Detection Efficiency 

The probability rj(p, 9 c . m .) that an elastic pp-scattering event 
leads to a valid trigger and survives all offline cuts was de- 
termined by Monte-Carlo simulations. It uses the available in- 
formation on the vertex-distribution, beam position and angles 
and accounts for the electronic thresholds and trigger condi- 
tions. All simulated events are then analyzed by the same soft- 
ware as used for the experimental data, applying the same re- 
construction algorithms and cuts. The fraction of events pass- 
ing these cuts r)(p, c .m.)> is fairly constant between 0.94 and 

0. 96 and shows no dramatic momentum dependence. The ef- 
ficiency is smallest for the lowest momenta and forward scat- 
tering and for the highest momenta and symmetric scattering, 

1. e. #iab.i ~ #iab,2 or c . m . « 90° . The main reduction is 
due to secondary reactions of the ejectiles in the beam-pipe or 
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Fig. 14. Comparison of two excitation functions with (closed) and without (open symbols) using the fuzzy-filter approach for background 
subtraction. To remove the main energy dependence all cross sections have been divided by the predictions of the PSA-solution FA00 of 1 8]. 



the detector, so that the information on one or both protons is 
sufficiently distorted to be removed on the trigger-level or by 
offline-cuts. The loss of 4 to 6 % is consistent with the expecta- 
tion derived from the known total hadronic cross-sections and 
the thickness and composition of matter in the detector set-up. 
For application, rj(p,9 CAn ) is approximated by a polynomial 
expansion in cos# cm . with momentum dependent coefficients. 

3.8 Detector Acceptance 

To account for the acceptance limits of the EDDA detector two 
software cuts, on 6> c . m . and </>, are applied in the analysis: The 
limits of the acceptance in the cm. scattering angle c . m . is 
given by the requirement that both protons are detected in the 
half-ring layer of the EDDA-detector. The minimum scattering 
angle accepted for the final data (shown as the solid line in 
Fig.GTJ was selected large enough to rule out elastic events to 
be lost due to the extended beam-target overlap or small-angle 
scattering of the ejectiles in the beampipe and the detector. 

Although the EDDA detector yields full azimuthal cover- 
age a cut on <fi is applied as well: For <\> close to 90° or 270°, 
events are lost on the trigger level when both particles hit the 
same - left or right - half of the ring layer. In addition the 
angular resolution achieved with the semi-rings is reduced at 
these angles, since close to the readout large variations in the 
light collection efficiency occur. Thus, we artificially reduced 
the azimuthal acceptance by selecting only events with the left 
particle emitted at |^>| < 76° and accounted for this reduced 
acceptance with a factor 180/152 = 1.184±0.007in the cross- 
section calculation. 

3.9 Normalization 

In internal target experiments the luminosity depends on the 
overlap of the target density distribution with the transverse 



profile of the stored beam. During acceleration both the hori- 
zontal beam position and the emittance - determining the beam 
width - vary considerably (cf. Fig. |6) and cannot be deter- 
mined, e.g. from elastic scattering data (cf. Sect. 13.51 . with 
sufficient precision. Therefore, two electromagnetic processes 
with known dependence on the instantaneous luminosity and 
beam energy were recorded concurrently with the elastic scat- 
tering data. These are the current of secondary electrons (SEM) 
emanating from the fiber target which is a function of the en- 
ergy deposit by the beam, and the rate of elastically scattered 
electrons, so-called ^-electrons, in PIN-diodes placed at 40° 
behind thin windows in the beam pipe (cf. Fig.[3J. These two 
methods are referred to as SEM- and PIN-monitors. 

As it will be outlined in the following sections, both are 
not accurate enough to allow an absolute luminosity measure- 
ment. However, the main uncertainties enter as multiplicative 
factors which are independent of beam-momentum and can be 
eliminated by an absolute normalization to high-precision ref- 
erence data 1 30] at one beam momentum (p rc f = 1 .455 GeV/c). 
The change of the luminosity with respect to this reference mo- 
mentum is provided by the SEM- and PIN-monitors with high 
accuracy, and will be referred to as relative luminosity deter- 
mination. 

Furthermore, both monitors are sensitive to the electron 
contents of the target only and measure L e , i.e. the luminos- 
ity with respect to proton-electron interactions. To relate this to 
the Luminosity Lh for proton-proton interactions the hydrogen 
~p and electron p e densities in the target enter, viz 

r p(p) T P , x Po T n ~ 

Lh = L e = — (jp) L e , (15) 

Pe Po 

Lsem or LpiN 

where po is the hydrogen density of an unused CH2-target. 
The hydrogen density ~p is the weighted average over the re- 
gion of the target sampled by the COSY-beam and therefore 
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momentum dependent. Its value will gradually decrease due to 
radiation damage in the course of the experiment and will be 
discussed in Sect. 13.9.31 For convenience the factor po/p e has 
been included in £sem and -Lpin, so that it is fixed by the abso- 
lute normalization and the factor — (»), common to both moni- 
tors, is deduced separately. A corresponding, but much smaller 
change in the electron density will also be taken into account. 



3.9.1 Secondary Electron Monitor (SEM) 

The SEM rate is proportional to the average energy AE de- 
posited in the target fiber ll3~ll . and therfore to the restricted 



energy loss rate 



dE 
dz 



T<T C , 



that energetic knock-on electrons above some energy T cut will 
escape from the target and not contribute to the energy deposit. 
We fix the cut-off energy to 9 keV by requiring that the effec- 
tive range of electrons equals half the target thickness. How- 
ever, any choice within a factor of 5 yields the same energy 
dependence within 0.6%, but not the same absolute value of 
the restricted energy loss. Thus, we expect the current Isem 
of emitted secondary electrons, as recorded by a sensitive am- 
peremeter [21], to be given by 



-?SEM = fcsEM • -^SEm(p) 



dE 






T<T cut 



(16) 



The current is carried by the electrons emitted from the target 
surface, and is dominated by electrons of very low energy, pro- 
duced by thermalization of the deposited energy. Hence, the 
proportionality factor fcsEM should not depend on the precise 
initial distribution of the electron energies and thus the beam 
momentum - it will, however, be sensitive to the structure of 
the surface. Any local damage of the aluminum-coating of the 
fiber and a resulting possible charge-buildup would make fcsEM 
a function of the position x along the fiber. Since the proton 
beam moves along the target horizontally during acceleration, 
this would result in an indirect dependence of fcsEM on p. Such 
an effect has been ruled out by comparing -Lsem to the elas- 
tic pp-scattering rate when steering the beam across the target 
horizontally during the flattop (cf. Fig.[6j, so that p is fixed and 
only x was varied. 



3.9.2 5-Electron (PIN) Monitor 

Elastically scattered electrons are detected at 40° by 500 /im 
thick PIN diodes with an area of 10 x 10 mm 2 . They are mounted 
in pockets at <\> = 0° and 180° in the beampipe, behind thin 
(250 /im) aluminum windows. The electrons are produced ac- 
cording to the well-known Rosenbluth cross section 1331 13411 
with the simplification, that due to the small momentum trans- 
fers («1 MeV/c) at COSY-energies the form-factors are unity. 
Detection at 40° offers both statistical precision ( da/df2\ a b ~ 
200 mb/sr) and electrons of sufficiently large energies (0.5- 
1 .5 MeV) to be observed with high efficiency. 

The luminosity is determined in two steps: experimentally 
the rate of scattered electrons N e is deduced from the singles 
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rate in the PIN-diodes Npw, corrected for the rate of hadrons 
N h . 

N e = N P1N - N h (17) 

Since all electrons are stopped within the PIN-diodes or their 
mounting, the rate of hadrons is determined by looking for co- 
incident hits in both bars and semi-rings of the outer detector 
layer (cf. Fig. [0 positioned behind the respective PIN-diode 
when viewed from the target. We found that 4-6% of all PIN- 
triggers had to be attributed to hadrons, with the fraction lin- 
early increasing with beam momentum. To estimate the number 
of hadrons escaping detection in the outer layer because of sec- 
ondary reactions or stopping, a telescope was mounted in place 
of a PIN-diode in an independent measurement. The telescope 
comprised a 5x5 mm 2 PIN-diode followed by a Cu-absorber, 
thick enough (0.7 mm) to stop all electrons, and a 10 x 10 mm 2 
PIN-diode to detect the hadrons. It could be shown that that 
0.5% of all PIN triggers not vetoed by the outer-detector lay- 
ers are due to hadrons |35] and that this fraction is not beam- 
momentum dependent. 

Secondly, the number of electrons has to be related by a 
Monte-Carlo simulation to the luminosity (cf. Eq. (IT5l ) 



Lpm(p) = —N e {p) ■ l M o(p) 

Pe 



(18) 



where po/Pe is the ratio of free protons to electrons in the tar- 
get and Imc is the average luminosity in p-e elastic scattering 
needed to produce one detected electron. The latter is given by 
MC-integration of the Rosenbluth cross-section, taking into ac- 
count the exact geometry of the PIN-diodes, all materials used 
for mounting and shielding, the COSY-beam direction, loca- 
tion and width (Sect. 13. 5> and the detection thresholds deter- 
mined experimentally. The interaction of the the electrons were 
modeled using the electron-gamma shower code EGS4 1 36J. In 
Fig-Elthe shape of the the experimental energy-loss spectrum 
is compared to MC-simulations, showing the transition from 
stopped electrons at lower energies to a typical Landau-like 
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Fig. 16. Relative change of the hydrogen contents of CH2-targets as 
a function of the accumulated dose as determined from measurements 
at fixed energies. The solid line is a parameterization as outlined in the 
text (from 1331 1 



energy-loss distribution at higher energies. Note, that - apart 
from normalizing to the same number of counts - no param- 
eters have been adjusted using the experimental data. The ab- 
solute energy scale was fixed by calibration with conversion- 
electrons from a 175 Hf-source. The main source of systematic 
errors in (mc 13711 are uncertainties in 

1 . the target location along the beam, 

2. the horizontal beam position, 

3. the thickness of the aluminum window 

4. and the precise size of the depleted area of the PIN-diodes. 

By taking the arithmetic mean of luminosity determined with 
the PIN-diodes mounted to the left and right of the beam, the 
second errors cancel and the last errors do not contribute when 
looking only at changes of ^mc with beam momentum. The re- 
sulting systematic error of Lpjn(p) / £piN(?>ref) increases with 
\p — Pref | and stays below 2.9% (1%) for beam momenta below 
(above) p Ici = 1.455 GeV/c. 

In principle an absolute luminosity could be deduced, but 
its precision hinges on the knowledge of the ratio of electrons 
to protons in the CH2-target. The naive estimate of po / p e be- 
ing four is changed by the aluminum coating (w20 pg/cm 2 ) 
to about 4.7. In addition, radiation damage (see Sect. I3.9.3> 
reduces the number of hydrogen in the target, as outlined in 
the next section. The absolute normalization to reference data 
agrees within 5-10% with our estimates of ~p/p e , consistent 
with the errors due to uncertainties of the aluminum thickness, 
the size of the depleted area of the PIN diodes and the accumu- 
lated dose of a specific target. 



3.9.3 Correction for Radiation Damage of the Target 

The hydrogen density of the CH2 -target decreases during the 
course of the experiment. This is due to radiation-induced cross- 
linking of the polypropylene polymers by replacing two C-H 
bonds by a C-C bond and a H2 molecule emanating from the 
target. In comparison, losses due to hadronic interactions are 
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Fig. 17. Average relative hydrogen density as a function of momen- 
tum for three samples corresponding to different setups of the COSY 
beam. 



negligible. The relative change of the hydrogen density p, and 
to a lesser extent of the electron density, with respect to an un- 
damaged target po will be a function of the acquired dose. The 
dose is not deposited uniformly along the target, but reflects the 
density distribution of the beam over the course of many COSY 
machine cycles (cf. Fig.|6j. Finally, during the acceleration tar- 
get regions with different hydrogen density will be sampled at 
different momenta which needs to be corrected for. 

To this end the history of each target with respect to the ac- 
quired dose has to be closely monitored. The rate at which dose 
is acquired is related to the instantaneous luminosity Ln(t), 
the specific energy loss of beam particles in the target and the 
normalized horizontal beam profile P(x, t) viz 



AD 
~~&t 



dE 



(P(t)) 

T<T cut 



kL H (t)P(x,t) 
p(t)/po 



(19) 



where the constant k — 113.6 Gy mm mbcm/MeV encom- 
passes all information on target material and dimensions (see 
1331 for details). Since Lh is the luminosity with reference to 
the hydrogen content of the target, the factor p/po corrects for 
the average relative hydrogen contents sampled by the beam. 

The relative change of the hydrogen content has been deter- 
mined experimentally, by observing the relative change in the 
yield from pp-elastic scattering with respect to that of inelastic 
p-carbon reactions at constant momentum as a function of the 
acquired dose, calculated from the measured luminosity 1 38 ] . 
The result is shown in Fig. Altogether with a parameterization 
of the form 



Po 



-(D) =/exp(-Ai£>) + (l-/)exp(-A 2 D) 



(20) 



The loss is described by a fast process with decay constant 
Ai = (4.4 ± 1.3) ■ 10~ 8 Gy -1 which reduces the hydrogen 
density by at most / = (5.8 ± 0.7) %, and a slow process 
with decay constant A 2 = (5.67 ± 0.28) • lO" 10 Gy -1 , with an 
additional scale uncertainty of A 1,2 of 1 1%. 

For each target the luminosity was recorded for the time of 
exposure to the COSY-beam and the beam profile P(x,t) is 
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Fig. 18. Comparison of the differential cross section of the three groups of samples (distinguished by the symbols O, □, and A). Only statistical 
errors are shown. To remove the main energy dependence all cross sections have been divided by the PSA-prediction of 1 8] (solution FAOO). 



known from the reconstructed vertex distribution (cf. Sect. 13. 5> 
and the dose can be calculated with Eq. dl9t in a first approxi- 
mation by assuming p/po = 1 to obtain the dose distribution. 
Using the hydrogen density as given by Eq. J20b the dose cal- 
culation can be iterated until self-consistency is reached. 

Finally, the hydrogen density as a function of the beam mo- 
mentum during the acceleration, averaged over the beam pro- 
file and the total time of the measurement is calculated for each 
data sample. Three examples, corresponding to different setups 
of the COSY-beam are shown. It turned out that for about one 
third of the data this correction is not momentum dependent 
(dashed line in Fig.ll7l>. because the change of beam position 
was small compared to the width of the beam. For other sam- 
ples the correction amounts to typically 5-15%. 

The shape of the correction factor only depends on the setup 
of the COSY beam. The data samples have been taken with 
three different settings, each one comprising roughly one third 
of the data. In Fig. ED excitation functions obtained with these 
three beam setups are compared at two cm. scattering angles, 
showing perfect agreement within statistics. When these cross 
sections are integrated over the full detector acceptance of the 
detector, they agree much better than the 2.5% error corridor, 
as given by the relative luminosity determination. 

Note, that the PIN-monitor is affected indirectly by changes 
in electron density p e of the target (cf. Eq. ( I18H . However, most 
of the electrons are attached to the carbon and aluminum nu- 
clei and the applied correction is smaller by a factor of about 
4.8. The SEM monitor is not affected by cross linking, since 
the SEM-yield is determined by the electron density in the alu- 
minum coating of the target. 

3.9.4 Absolute Normalization 

The SEM and PIN monitors record the change of the luminos- 
ity with beam momentum very accurately. The absolute value 



is obtained by normalizing the pp-elastic scattering angular dis- 
tribution at one projectile momentum bin p rc f to precise refer- 
ence data of differential cross sections. We use the precision 
measurement at p ref = 1.455 GeV/c (T rcf = 793 MeV) by Si- 
mon et al. 13011 . where the absolute normalization uncertainty 
is given as < 1 %. By numerically integrating the data points 
of 1 30] we arrive at a reference cross section 

r 88 - 51 ° da 
<7 re f = 27T / d6» c . m . sin6> c . m .— —(0cm.) 

J 39.84° rtiic.m. (21) 

= 11.16 ± 0.02 stat - ± 0.11 syst - mb. 

to be used for the final absolute normalization. By numerically 
integrating the angular distribution at p rc f of this experiment 
over the same angular range, the values of &sem an ^ Po I Pe 
(Eqs. d!6l > and ( I18H are adjusted to yield the same integrated 
elastic scattering cross section er rc f . 

The absolute normalization receives an additional statisti- 
cal error from the statistical precision of the angular distribu- 
tion of our data samples at the normalization momentum p rc f . 
The procedure to minimize the influence of this error on the 
final result when combining all data samples is described in 
Sect. lTHl 

3.9.5 Comparison of Luminosity Monitors 

In Fig. ^] the result for the integrated luminosity from the 
SEM- and PIN-monitors are compared for one typical data sam- 
ple. The agreement is typically much better than 2.5%. The er- 
ror on the luminosity derived from e-p elastic scattering (PIN- 
monitor) is dominated by systematic uncertainties from the MC- 
simulations. Furthermore, below 1.1 GeV/c the (5-electrons at 
40° are so low in energy, that they can no longer be detected 
by the PIN-monitor. Therefore the luminosity obtained from 
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Fig. 19. Time-integrated luminosity for one out of 17 data samples 
from the two luminosity monitors. The dashed line shows a ±2.5% 
band around Lsem to test the consistency to Lpin. The absolute nor- 
malization scale has been fixed atp rc f = 1455 MeV/c. Note the sup- 
pressed zero on the vertical axis. 
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Fig. 20. Comparison of the SEM and PIN luminosity monitors aver- 
aged over all data samples. The errors are dominated by the uncer- 
tainty in l M c- 



the SEM-monitor is used in the analysis and the PIN-monitor 
is taken as a consistency check and to derive an error estimate 
on the relative luminosity determination. When averaged over 
all data samples the two luminosity monitors give consistent 
results (cf. Fig. I20i and deviations stay below 1.25% within 
known statistical and systematic errors. Since some samples 
show deviations as large as 2.5% we use this value as an error 
estimate of the systematic relative luminosity determination. 
This corridor is shown in Figs.ll9landl20las dotted lines. 

The total absolute luminosity summed over all data sam- 
ples is shown in Fig. The first data taking period covered 
the momentum range 1.1-3.4 GeV/c and the second one 0.7- 
3.3 GeV/c. The spikes close to the maximum energy are due 
to the reduced ramping speed in the transition to the flattop 
of the accelerator cycle. Note that the luminosity is more or 
less constant over the whole momentum range (cf. Figs.ll9land 
12 It . This shows that the loss of beam particles due to the fiber- 
target is compensated by the increase in beam current due to the 
rising revolution frequency during acceleration. This is in line 
with the finding 1 2 21 13911. that the heating of the stored beam by 
the fiber-target is counterbalanced by the adiabatic damping of 
the betatron oscillations evinced as a near-constant beam-width 
6x v (cf. Fig. l3.5> . 
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Fig. 21. Total integrated luminosity for all data. 



3.10 Beam Momentum Determination 

The beam momentum for each scattering event is determined 
from the relative time t of the event with respect to the start of 
the COSY-cycle. This allows to calculate the nominal momen- 
tum p(t) of the COSY-beam based on the mathematical model 
used for programming the function generators of all compo- 
nents of the COSY ring, especially the cavity's RF and the 
dipole-currents. 

The beam momentum only depends on the revolution fre- 
quency v and the closed-orbit length C of the stored beam: 



p = m p /3~fc = m p c 



(22) 



Since COSY operates on the first harmonic the revolution fre- 
quency is equal to j/rp- The RF-frequency was recorded dur- 
ing data taking, and perfect agreement with its nominal value 
was found. Time stamps were recorded with the data every 
2.5 ms using a high-precision 20 MHz clock (relative accuracy 
of 10~ 6 ). Thus, based on the ramping speed of 1.15 MeV/c 
per ms, the momentum for individual events is known bet- 
ter than ±1.5 MeV/c. The remaining uncertainty is due to 
the closed-orbit length C and its possible momentum depen- 
dence. Upper limits on its deviation from the ideal orbit length 
(Ci = 183.472 m) were deduced from measurements of the 
beam positions as a function of time at 29 locations around 
the ring. With the parameters of the COSY-lattice the maximal 
change in C compatible with constraints from ion-optics was 
deduced. The corresponding relative momentum deviation is at 
most 3 • 10" 4 , i.e. less than 1 MeV/c at 3.4 GeV/c A 
correction of C with respect to C\ turned out to be negligible. 

3.11 Dead-Time Correction 

The dead time fraction of the data-acquisition system was typ- 
ically 90% and its instantaneous value had to be determined 
accurately. Using fast (20 MHz) scalers we recorded both the 
number of events read-out Nr and the number of all trigger Nt 
(about 10 kHz) dead-time free. The dead time fraction is then 
given by r(t) = (1 — Nr/Nt) with statistical uncertainties of 
less than 0.5 %. Systematic errors were checked by taking data 
samples with data rates varying by orders of magnitude and 
turned out to be negligible. 
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3.12 Combination of Data Samples and Consistency 
Checks 

All 17 data samples were analyzed separately, the relative lu- 
minosity was fixed using the SEM-monitor and the absolute 
cross section normalization with respect to 1 30] as described in 
Sect. l3.9.4l These data samples are distinguished by at least one 
of the following: the time of the measurement, the CH2-target 
used, COSY-beam parameters and trigger conditions. This al- 
lows to test for systematic variations of the results with these 
parameters. Let us introduce the shorthand-notations er,; and 
8o~i for the differential cross section of the i'th data sample and 
its error. The compatibility of the results {o~i} of n data sam- 
ples, i.e. i = 1, . . . , n, was tested by looking at the individual 
contributions to the total \ 2 when minimizing 



X 



EE 

»=1 .7>fe 



(M (Ti(9j,pk) - <Jj,k) 2 
M 2 6* 2 (9 3 , Pk ) 



(23) 



by variation of the normalization factors Mi with Mi fixed to 
unity. Here, the mean cross section is then given by the weighted 
mean 



E 



Mi (Ti(6j,p k 




1 



^ M 2 5a 2 {9 3 , Pk ) l ^ NfSa^p*) 



(24) 



The Mi were treated as free parameters, so that the statistical 
accuracy 8cr 2 (9j,p rc {) at the normalization momentum, dom- 
inating the absolute normalization error of the individual data 
samples, do not contribute to the \ 2 - The fitted values for the 
Mi turn out to be very close to unity, well within the uncertainty 
due to Sa 2 (p TC {). 

When all data samples were combined the total \ 2 P er de- 
gree of freedom was 1 .03. When testing the influence of certain 
aspects of the analysis on the final results, samples with the 
same conditions with respect to the aspect under study were 
combined viz Eq. 1231 and renormalized to l30l in order to 
increase the statistical accuracy of the consistency check. An 
example is shown Fig.[D3]where the dependence on the COSY- 
beam setup is tested, showing consistent results within the sta- 
tistical errors. No systematic deviation of the results from dif- 
ferent samples could be determined. In addition cuts applied to 
the data have been varied within reasonable limits to check for 
possible systematic dependencies of the final result on details 
of the analysis. As an example, the value of a cut , the maximum 
allowed value for the kinematic deficit may be increased by 3° 
without changing the results significantly, although at the ex- 
pense of increased statistical error. The change is always much 
smaller than the estimated systematic error and amounts to less 
than 1% for the majority of the data. 

Final cross sections are obtained by combining all 17 sam- 
ples with the help of Eq. ( 1231 and renormalizing them at p Te { to 
1 30]. These data comprise a total of 37- 10 6 pp elastic scattering 
events. 



3.13 Errors and Uncertainties 

The calculation of the excitation functions receives error con- 
tributions from all factors on the right hand side of Eq. l|3}. We 
distinguish five different kinds of errors: 



1 . statistical errors, due to the number of observed counts from 
CH2- and C-targets as in Eq. l|8) with a very small contri- 
bution from the relative luminosity L of the carbon-sample. 
This error ranges from 1% at small 6* cm . and small mo- 
menta to 7% at 9 C . m .~ 90° and large momenta reflecting 
the functional dependence of dcr./df2 c . m . 

2. systematic errors which are treated as uncorrected between 
data points adjacent in scattering angle or momentum. Here, 
artefacts from the reconstruction of the scattering angles 
contribute. The angular resolution of the EDDA detector is 
not entirely homogeneous, but varies slightly with Q\i (cf. 
Eq. 0). When the data are binned in 9 c . m , an uncertainty of 
2.4%, usually anti-correlated for neighboring bins, arises. 
Smaller contributions to this uncertainty stems from the cut 
on the azimuthal angle tj> (0.6%) and the loss of events due 
to the cut on P oL (0.5%). 

3. systematic errors associated to undetected elastic scatter- 
ing events and misidentified inelastic events. This uncer- 
tainty will have a smooth dependence on scattering angle 
and beam momentum. The detection efficiency, obtained 
by MC methods, contributes with 1.5% while possible con- 
tamination with inelastic events account for an uncertainty 
increasing with beam momentum and scattering angle from 
less than 1% to 5.5%. These contributions are always smaller 
than the first two errors combined. 

4. a relative normalization error, common to all data points at 
the same momentum of 2.5%, derived from the maximum 
discrepancy observed in the relative luminosity determina- 
tion of the SEM and PIN monitors. The statistical error of 
the SEM signal as well the uncertainty in determination of 
the DAQ dead-time are negligible in comparison. Uncer- 
tainties in the correction ~p/po, due to uncertainties in A1.2 
of Eq. J20i and the applied dose, enter only as the ratio to 
its value at p re f, so that its contribution is less than 0.3% for 
all beam momenta. 

5. an absolute normalization error common to all data of less 
than 1.5%, comprising the uncertainty of <r re f of Eq. 11211 
and the statistical errors of our data at p re f . 

In all figures presenting the data of this work, only the first two 
errors are shown, i.e. the statistical error and an additional 2.5% 
systematic error. All errors are listed in B41I1 . 



4 Experimental Results 

The experimental results of the present elastic pp scattering ex- 
periment are the unpolarized differential cross sections der / dil c l 
as a function of the c.m.-scattering angle 9 c . m . and the labora- 
tory momentum p of the proton beam. Each differential cross 
section refers to bins A6 c . m . = 2° and Ap = 25 MeV/c, 
which were chosen to have reasonable statistical precision for 
each data point. We covered a beam momentum range of 0.7 — 
3.4 GeV/c corresponding to laboratory kinetic energies 0.23 — 
2.59 GeV and total cm. energies ^/s = 2.0 — 2.9 GeV. Since in 
pp elastic scattering the final state particles are indistinguish- 
able and cross sections are symmetric with respect to 9 c . m , = 
90°, data are only given for the cm. angular range 34° — 90°. 
All 2888 data points are available online through the world- 
wide web 1 41]. They may be viewed as 28 excitation functions 
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Fig. 22. Comparison of the differential cross section of 1 14] (open symbols) with this work (closed symbols) at two angles. To remove the main 
energy dependence all cross sections have been divided by the PSA-prediction of 1 8] (solution FAOO). 



(cf. Fig. I23> at cm. scattering angles between 34° and 90° 
or 108 angular distributions (cf. Fig.l25> at 108 beam momenta 
between 0.7 and 3.4 GeV/c. Note, that due to kinematics the an- 
gular acceptance is reduced at very small and very high beam 
momenta. 

As compared to our previously published results 1 14], based 
on a subset of the data, the momentum range is increased from 
1.1-3.3 GeV/c to 0.7-3.4 GeV/c and the statistical precision in- 
creased by up to a factor of 3. In addition, in the analysis of 
Ref. [14] the variation of the hydrogen density along the tar- 
get due to radiation damage was not corrected for, this led to 
a systematic underestimation of the cross section by about 5% 
above 2 GeV/c, as displayed in Fig. [22] Note, that the results 
from 1 14] are superseded by this work and should not be used 
anymore. 



4.1 Comparison to Other Data 

In Fig. |24| the available database on pp elastic scattering (cf. 
I7ll8ll42ll43ll44ll and references therein) is plotted for the same 
angles as for the data of the present work in Fig.|23] The ben- 
efit of a consistent normalization, made possible by measuring 
during acceleration in an internal experiment, as well as the 
improved statistical accuracy is evident. Previous data scatter 
considerably around available phase-shift solutions (discussed 
in detail in Sect. I5.2> partly due to larger statistical errors, but 
mainly due to differences in the absolute cross section normal- 
ization of the various experiments. Most notably the data of 
Jenkins et al. I45II from the ZGS are lower by about 20% (irA 
in Fig.|24> and clearly disagree with our data as well as those 
of 1 46l I47I1 . In Fig. [26] angular distributions are compared at 
four momenta to data from other experiments. The shapes of 
the differential cross section are consistent within quoted un- 
certainties, however, again absolute normalizations are at vari- 
ance. To show the size of this scale difference we have renor- 



malized our data using the procedure described in Sect. l3.9.4l to 
match the normalization of other experiments which spanned 
a larger range of beam momenta. The result is displayed in 
Fig. |^]and shows that the data of Kammerud et al. I46II are 
consistently higher by 10 %, the data of Williams et al. |47] 
agree, within the sizeable uncertainties, and the data from Al- 
brow et al. |48] are partly compatible with our normalization. 
On the average these data tend to be larger at momenta above 
about 1 .5 GeV/c, with the exception of the data of Jenkins et al. 
1 45] discussed earlier. The excitation function of Garcon et al. 
H49I1 at 9 C . m . « 90° , measured by a similar technique, shows 
the opposite trend (Fig. [28} and cross sections are smaller by 
up to 2 standard deviations. 

Note, that all 2888 cross section data obtained in this exper- 
iment share a common absolute normalization factor. In case 
data becomes available in the future, with a superior method to 
obtain an absolute normalization, our data could be renormal- 
ized by a common factor to be determined along the same lines 
as described in Sect. 13.9.41 The relative normalization uncer- 
tainty of 2.5 % would not be affected. 



5 Discussion 

Theoretical models of the nucleon-nucleon interaction have been 
very successful in describing elastic scattering data at energies 
where cross sections of inelastic reactions are still small. In this 
domain (T p below about 0.5 GeV) high-precision potentials, 
purely phenomenological or based on the meson-exchange pic- 
ture, are available (for a review see Ref. |62]). More recently 
effective field-theory has entered the stage and produced re- 
sults of comparable quality, cf. Refs. |63, 64, 65] and references 
therein. However, in the energy region between 0.5 GeV and a 
few GeV little progress has been made in recent years. In the 
late 80s, the meson exchange models were extended to higher 
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Fig. 23. Excitation functions of do-/df2 c . m . for six cm. scattering an- 
gles in comparison to SAID solutions SM94 (dotted), SM97 (dashed) 
and FA00 (solid line). 



Fig. 24. Collection of published data (From the database of 1 8] without 
EDDA data of Albers et al. [ 14]) plotted as excitation functions at six 
angles (±1°) in comparison to phase shift predictions as in Fie. 1231 



energies by including inelastic reactions by coupling to inter- 
mediate NA and A A channels i6rll67ll . A qualitative descrip- 
tion of the NN-data below 1 GeV has been achieved. However, 
at higher energies these models fail badly and differential cross 
sections at large scattering angles are strongly overestimated 
ir38il . This may not be surprising since inelasticities at these en- 
ergies will no longer be driven predominantly by the Z\(1232) 
resonance. Work aiming to improve these models are currently 
underway Iffijl . 

On the high-energy side optical potential models and Regge- 
phenomenology has been used to describe the gross features of 
the unpolarized differential cross section 14 51 14(11. Dimensional 
scaling 1 70] predicts cross sections to follow a s~ 10 behavior at 
large momentum transfer. Our data are too low in momentum 
transfer for these models to be strictly applicable. Nonetheless, 
in Fig. [29] we compare our differential cross sections, plotted 
as ^n-vs — t, to the prediction of 1 70] (solid line, note that for 



dt 

9c.m.= 90° 



we have s 



Ami 



2t), showing that we may 



be barely touching the region of dimensional scaling, in agree- 
ment with 1 45 ] . 

The immediate benefit of the high-precision data obtained 
in this experiment is two- fold: first, its consistent normalization 
allow for stringent tests of energy-dependent structures, as they 
could arise from coupling to an intermediate resonant state, and 
secondly it will further help to consolidate phase-shift parame- 
ters for the isotriplet elastic NN-channel. 



5.1 Upper Limits on Resonant Contributions 

All excitation functions show a smooth and rather structureless 
dependence on beam momentum. No sharp energy-dependent 
structure is observed which could be taken as evidence for a 
narrow resonance. Therefore, upper limits for the elasticities 
r e i /r of hypothetical narrow resonances are deduced from the 
smooth excitation functions. To this end a Breit-Wigner reso- 
nance term is introduced into the scattering matrix element of 
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Fig. 25. Angular distributions of dcr/dfic.m. (•: data of this work) for 
ten beam momenta p in in comparison to SAID solutions (cf. Fig. l23i . 



the partial wave which is assumed to exhibit resonance behav- 
ior. The interference between the Breit-Wigner resonance term 
and the non-resonant amplitudes determines the size of a reso- 
nance excursion in the excitation functions. The non-resonant 
amplitudes represent the null-hypothesis. In order to establish 
the null-hypothesis the EDDA data are fitted by a special en- 
ergy dependent phase shift analysis along Ref.|8|. 

In the test calculations the hypothetical resonance energy 
E R , total width r and resonance phase (f> R are varied system- 
atically in the range E R = 2.2 - 2.8 GeV, r =10— 100 MeV 
and <f) R = Q° — 360° . The hypothesis of the existence of a res- 
onance in a partial wave is tested by gradually increasing the 
partial elastic width F e i until the resonance is excluded within 
99% confidence level by a x 2 -test based on the the differential 
cross sections and analyzing power data Hal . For the unknown 
phase tpji the value giving the largest, and thus most conser- 
vative limit on F e i /T is chosen. Typical upper limits of F e i /T 
for the five lowest uncoupled partial waves are 0.08 for 1 Sq, 



Fig. 26. Angular distributions at four beam momenta. Data of this 
work (•) are compared to published data, where the symbols O, □, 
A, 0,O, and )K correspond to data from 1^ l5(i0ll R^. pjh at 
1.08 GeV/c, IM^I^lHlsllsl^ at 1.45 GeV/c, j43|43|5J at 
2 GeV/c and [ 43l^ll^alHalHlllEalH3] at 3 GeV/c. 



0.04 for ^2, 0.10 for 3 P , 0.03 for 3 Pi and 0.05 for 3 F 3 . For 
instance the 1 So dibaryon resonance predicted by Lomon et al. 
IToft at E R = 2.7 GeV with r el /T = 0.1 and 7" = 50 MeV 
can safely be excluded. The method is described in detail in a 
forthcoming paper 1 26] . 

It should be noted in this context that two broad resonant 
structures are well known at lower energies, the 1 D2 resonance 
at E R w 2.15 GeV and the 3 F 3 resonance at E R » 2.17 GeV 
with widths of around 120 MeV. They appear as counterclock- 
wise circles in the Argand diagrams of the phase shift analysis 
lE2ll of the world data set. Especially the excitation functions 
of the spin-dependent total cross sections Aul and Aar mea- 
sured at ZGS I7ll l72ll show marked energy-dependent struc- 
tures. These X D2 and 3 F3 resonances in the pp-system can be 
interpreted conventionally by intermediate 5 S2 and 5 P3 NA 



D. Albers et al.: A Precision Measurement of pp Elastic Scattering Cross Sections at Intermediate Energies 



21 





1.3 




1.2 






C 


1.1 


O 




U 


1 










0.9 


. N 




s 


0.8 


c 
o 


0.7 


c 






0.6 




a Williams '72 

* Jenkins '80 

• Simon '93 



1.5 2 2.5 
p (GeV/c) 



Fig. 27. Comparison of the absolute normalization constant relative 
to |3(J at 1.455 GeV/c for angular distributions measured at CERN 
(Albrow et. al [48]), the Argonne ZGS (Kammerud et. al 14^11. R ufher- 
ford (Williams et al. |47]), Brookhaven (Jenkins et. al 1 45]), and 
LAMPF (Simon et al. 13(11 ). The smaller error bars show statistical 
uncertainties only and the larger error bars contain the normalization 
uncertainty added in quadrature. The dashed horizontal line shows the 
relative normalization error of this experiment. 
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Fig. 28. Comparison to the excitation function of 1 49] at # c . m . ~ 90° . 
The EDDA date have been interpolated to the exact angle as in |49] 
and all cross sections have been divided by d(jfl t /df2 c . m ., a smooth 
parameterization of the data. 



states These broad resonances show up as 

a rather rapid change in the shape of the angular distribution 
around 1.25 GeV/c {y/s = 2.15 GeV) where the cross section 
exhibit a steep ascend (descend) at forward (backward) angles. 



5.2 Impact on phase shift analysis 

The excitation functions of the differential cross sections are 
measured over a wide momentum range from 0.7 to 3.4 GeV/c. 
They have a great impact on the energy-dependent phase shift 
analysis since they represent a precise and consistent data set 
of 2888 data points. Another important fact is that the unpolar- 



t (GeV ) 



Fig. 29. Differential cross section at 3 selected cm. scattering an- 
gles: 89° (□), 69° (O), and 49° (0). 



ized differential cross sections fixes the global scale of all am- 
plitudes and thus enters repeatedly into the phase shift analysis 
through all spin dependent cross sections, which are products 
of spin-observables with the unpolarized cross section. 

The dotted curves in Fig.[23]are phase shift solutions SM94 
of Arndt et al. 1 6] from summer 1994, i.e. before the first EDDA 
data were available. The maximum kinetic energy of this phase 
shift solution was 1 .6 GeV corresponding to a beam momen- 
tum of 2.36 GeV/c. After the publication of the first EDDA data 
[14] the VPI group extended their energy dependent phase shift 
analysis from 1 .6 GeV (2.36 GeV/c) up to 2.5 GeV (3.3 GeV/c), 
with the solution SM97. Meanwhile an energy dependent phase 
shift solution FA00 is available up to 3.0 GeV (3.82 GeV/c) 
laboratory kinetic energy (beam momentum) The phase 
shift solutions SM97 and FA00 are shown as dashed and solid 
curves in Figs. |23] and |25], respectively. Phase-shift predictions 
show an apparent oscillation both in angle and beam momen- 
tum about the data of this work. This may be an artefact of the 
current parameterization of the energy-dependence of phase- 
shifts. Therefore, including the data of the present work will 
not only slighly modify the phase-shifts to reflect the differ- 
ences to our previous results ( II14II and Fig. I22> but may also 
allow to improve the ansatz for the variation of the phase shift 
parameters with energy. In this process a better description of 
the angular distribution at the higher momenta should be at- 
tempted. 



6 Summary 

In a dedicated experiment, protons accelerated in a synchrotron 
have been scattered of internal CH2 -fiber targets during acce- 
laration at beam momenta between 0.7 and 3.4 GeV/c (T p = 
0.23 . . . 2.59 GeV). Elastic scattering events have been identi- 
fied by the EDDA-detector over a wide angular range (# c . m . = 
34° . . . 90° ). This experimental technique allows a precise 
monitoring of the relative change in luminosity with beam mo- 
mentum which leads to a consistent normalization of data taken 
at different momenta and is therefore ideally suited to measure 
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excitation functions with high precision. The absolute cross 
section scale was fixed at the reference momentum 1 .455 Ge V/c 
to high-precision data from LAMPF 13011 . The average com- 
bined statistical and systematic error is 3.8%, it ranges from 
3% for small momenta and scattering angles to 10% at high 
momenta and 6> c m . = 90° . The total absolute normalization 
uncertainty of the complete data set is below 1.5% with an ad- 
ditional relative normalization uncertainty of 2.5% common to 
all data points at the same beam momentum. 

The wealth of new data obtained in this experiments re- 
places the data reported earlier I14ll and utilized improved anal- 
ysis methods to yield higher statistical and systematic preci- 
sion, most notably at higher energies. These data will further 
improve extracted phase-shift parameters and provide an im- 
portant normalization standard for measurements of NN reac- 
tions. The data are available via the world-wide web [41]. 

The excitation-functions obtained with steps of 25 Me V/c 
in beam-momentum (or A*Js « 8.5 MeV) have been applied 
in tests for contributions of narrow resonances, like dibaryons, 
coupling to the elastic channel. No evidence for any narrow 
(10-100 MeV) structures have been found in the invariant mass 
range 2.0 ... 2.8 GeV covered by the experiment, ruling out 
a strong coupling of such resonances - if they exist - to the 
elastic channel. 

We thank the operating team of COSY for excellent beam support. 
This work was supported by the German BMBF, contracts 06BN664I(6) 
and 06HH852, and by the Forschungszentrum Julich (Auftrage 411 
268 03 and 41 126903). 



References 

1. C. Lechanoine-LeLuc and F. Lehar, Rev. Mod. Phys. 65 
(1993). 

2. W. Haeberli et al., Phys. Rev. C55, 597 (1997). 

3. B. von Przewoski et al., Phys. Rev. C58, 1897 (1998). 

4. F. Rathmann et al., Phys. Rev. C58, 658 (1998). 

5. J. Ball et al., CTU Reports 4, 3 (2000). 

6. R. A. Arndt, I. I. Strakovskii, and R. L. Workman, Phys. 
Rev. C50, 2731 (1994). 

7. R. A. Arndt, C. H. Oh, I. I. Strakovsky, R. L. Workman, 
and F. Dohrmann, Phys. Rev. C56, 3005 (1997), nucl- 
fh/9706003. 

8. R. A. Arndt, I. I. Strakovsky, and R. L. Workman, Phys. 
Rev. C 62, 34005 (2000). 

9. J. Ball et al., Phys. Lett. B 320, 206 (1994). 

10. P. Gonzalez, P. LaFrance, and E. L. Lomon, Phys. Rev. 
D35, 2142 (1987). 

11. P. J. G. Mulders, A. T. M. Aerts, and J. J. de Swart, Phys. 
Rev. Lett. 40, 1543 (1978). 

12. P. J. Mulders, A. T. Aerts, and J. J. D. Swart, Phys. Rev. 
D21, 2653 (1980). 

13. P. LaFrance and E. L. Lomon, Phys. Rev. D34, 1341 
(1986). 

14. D. Albers et al. (EDDA), Phys. Rev. Lett. 78, 1652 (1997). 

15. M. Altmeier et al. (EDDA), Phys. Rev. Lett. 85, 1819 
(2000). 

16. F. Bauer et al. (EDDA), Phys. Rev. Lett. 90, 142301 
(2003). 



17. R. Maier, Nucl. Instr. and Meth. A390, 1 (1997). 

18. K. Ackerstaff et al. (EDDA), Nucl. Instrum. Meth. A335, 
113 (1993). 

19. J. Bisplinghoff et al. (EDDA), Nucl. Instr. and Meth. A329, 

151 (1993). 

20. D. Albers et al. (EDDA), Nucl. Instr. and Meth. A 371, 388 
(1996). 

21. H. Scheid, Dissertation, Institut fur Strahlen- und Kern- 
physik, Universitat Bonn (1994), 

see: http ://w w w. iskp . uni-bonn.de/edda/dipldiss. html 

22. F. Hinterberger and D. Prasuhn, Nucl. Instr. and Meth. 
A279,413 (1989). 

23. J. Lindlein, Diplomarbeit, I. Exp. Phys., Universitat 

Ha mburg (1996), 

see, http://kaa.desy.de/edda/papers/Diplomalist.inhalt.html 

24. J. P. Coffin, Experimental Techniques in Nuclear Physics 
(de Gruyter, Berlin, 1997), p. 309. 

25. K. Ackerstaff et al. (EDDA), Nucl. Instrum. Meth. A491, 
492 (2002). 

26. H. RohdjeB et al. (EDDA) , nucl-ex/0403043. 

27. H.-H. Bothe, Fuzzy Logic - Einfiihrung in Theorie undAn- 
wendungen (Springer, Heidelberg, 1995). 

28. F. Hoppner, F. Klawonn, R. Kruse, and T. Runkler, Fuzzy 
Cluster Analysis (Wiley, New York, 1999). 

29. M. Busch, Dissertation, ISKP, Universitat Bonn (2001) 
see: http://www.iskp.uni-bonn.de/edda/dipldiss.html 

30. A. J. Simon et al., Phys. Rev. C 48, 662 (1993). 

31. E. J. Sternglass, Phys. Rev. 108, 1 (1957). 

32. Particle-Data-Group, C. Caso et al. (Particle Data Group), 
Eur. Phys. J. C3, 1 (1998), see also http://pdg.lbl.gov 

33. M. N. Rosenbluth, Phys. Rev. 79, 615 (1950). 

34. G. Kallen, Elementary Particle Physics (Addison- Wesley, 
1964). 

35. T. Huskes, Diplomarbeit, ISKP, Universitat Bonn (1997), 
see: http://www.iskp.uni-bonn.de/edda/dipldiss.html 

36. W. R. Nelson, H. Hirayama, and D. W. O. Rogers, The 
EGS4 Code System (1985), SLAC-Report 265, URL: 
' http://www.slac.stanford.edu/egs/ '. 

37. R. GroB-Hardt, Dissertation, Institut fur Strahlen- und 
Kernphysik, Universitat Bonn (2001). 

38. H. RohdjeB et al. (EDDA), nucl-ex/0403044. 

39. F. Hinterberger and D. Prasuhn, Nucl. Instr. and Meth. 
A321, 453 (1992). 

40. H. P. Engelhardt, Dissertation, Institut fur Strahlen- und 
Kernphysik, Universitat Bonn (1998), 

see: http ://w w w. iskp . uni-bonn.de/edda/dipldiss. html 

41. Data can be accessed via http://kaa.desy.de and 
http ://w w w. iskp . uni-bonn.de/gruppen/edda 

42. R. A. Arndt, L. D. Roper, R. L. Workman, and M. W. Mc- 
Naughton, Phys. Rev. D45, 3995 (1992). 

43. R. A. Arndt, J. M. Ford, and L. D. Roper, Phys. Rev. D32, 
1085 (1985). 

44. R. A. Arndt, W. J. Briscoe, R. L. Workman, and I. I. 
Strakovsky, SAID: Scattering Analysis Interactive Dialin 
(2004), http://gwdac.phys.gwu.edu/analysis/nnjnalysis. 

45. K. A. Jenkins et al., Phys. Rev. D 21, 2445 (1980). 

46. R. C. Kammerud et al., Phys. Rev. D 4, 1309 (1971). 

47. D. T. Williams et al., Nuovo Cimento 8A, 447 (1972). 

48. M. G. Albrow et al., Nucl. Phys. B 23, 445 (1970). 



D. Albers et al.: A Precision Measurement of pp Elastic Scattering Cross Sections at Intermediate Energies 



23 



49. M. Garcqn et al., Nucl. Phys. A 445, 669 (1985). 

50. K. Abe et al., Phys. Rev. D12, 1 (1975). 

51. A. J. Simon et al., Phys. Rev. C 53, 30 (1996). 

52. G. W. Hoffmann et al., Phys. Rev. C37, 397 (1988). 

53. D. Ottewell et al., Nucl. Phys. A412, 189 (1984). 

54. A. v. Dobrovolsky et al., Nucl. Phys. B214, 1 (1983). 

55. M. 1. Barlett et al., Phys. Rev. C27, 682 (1983). 

56. B. A. Ryan, A. Kanofsky, T. J. Devlin, R. E. Mischke, and 
P. F. Shepard, Phys. Rev. D3, 1 (1971). 

57. F. Shimizu et al., Nucl. Phys. A, 445 (1982). 

58. A. R. Clyde, Lawrence Livermore Radiation Laboratory 
Report No. UCRL-16275 (1966), unpublished. 

59. C. M. Ankenbrandt et al., Phys. Rev. 170, 1223 (1968). 

60. D. R. Rust et al., Phys. Rev. Lett. 24, 1361 (1970). 

61. I. Ambats et al., Phys. Rev. D 9, 1 179 (1974). 

62. R. Machleidt, in Adv. Nucl. Phys., edited by J. W. Negele 
and E. Vogt (Plenum Press, New York, 1989), vol. 19, 
chap. 2, pp. 189-376. 

63. U. van Kolck, Prog. Part. Nucl. Phys. 43, 337 (1999). 

64. E. Epelbaum, W. Glockle, and U.-G. Meissner, Nucl. Phys. 
A671, 295 (2000), nucl-th/99 10064. 

65. D. R. Entem and R. Machleidt, Phys. Rev. C68, 041001 
(2003),nucl-th/0304018. 

66. C. Elster, K. Holinde, D. Schiitte, and R. Machleidt, Phys. 
Rev. C38, 1828 (1988). 

67. C. Elster, W. Ferchlander, K. Holinde, D. Schiitte, and 
R. Machleidt, Phys. Rev. C37, 1647 (1988). 

68. K. O. Eyser, R. Machleidt, and W. Scobel (2003), nucl- 
th/03 11002. 

69. C. Elster, Ohio University, Athens, Ohio (2003), private 
communication. 

70. S. J. Brodsky and G. R. Farrar, Phys. Rev. Lett. 31, 1153 
(1973). 

71. LP. Auer et al., Phys. Rev. Lett. 41, 354 (1978). 

72. E. K. Biegert et al., Phys. Lett. 73B, 235 (1978). 

73. D. V. Bugg, J. Phys. G 10, 717 (1984). 

74. P. Gonzalez and E. L. Lomon, Phys. Rev. D34, 1351 
(1986). 

75. P. LaFrance, E. L. Lomon, and M. Aw (1993), nucl- 
th/9306026. 

76. T. Ueda, Phys. Lett. B119, 281 (1982). 



